火星表面热环境对航天器热控影响分析

火星大气对太阳辐射产生吸收和散射作用，同时还将与火星表面航天器发生对流换热。热设计时难以直接评估对流、辐射和导热三种换热对航天器的影响，从而确定主要的控温途径。在调研火星表面辐射、大气等热环境的基础上，从线性化传热系数和对流辐射比的角度对比分析了辐射、对流和导热对航天器的影响。器表辐射传热系数随光学属性和温度的变化范围为0.3～1.4W/(m2·℃)，对流传热系数随风速变化为0.2～1.5W/(m2·℃)，器内导热传热系数可控制在0.25W/(m2·℃)以下。结果表明，太阳辐射较火星表面和天空辐射而言是主要外热源，航天器表面的辐射和对流换热为两条并联换热途径，两者均可成为主要换热途径，器内导热传热是控制航天器内外隔热的主要可控因素。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

限带白噪声随机过程的雨流幅值概率密度函数模型

对17种不同带宽系数的限带白噪声随机过程进行了雨流循环计数统计,提出了一个限带白噪声随机过程的雨流幅值概率密度函数(Probability density function,PDF)模型。该模型是Rayleigh分布和Weibull分布的线性组合,其中待定系数均为随机过程谱参数的函数。采用该模型对这17种随机过程的雨流幅值概率密度函数进行了公式拟合,探究了模型中待定参数与随机过程谱参数之间的关系,确定了模型表达式。对照随机过程的雨流计数统计结果,将该模型与Dirlik模型的预测精度进行了比较,结果表明该模型的预测效果优于Dirlik模型。     

散射对高温气体-碳黑颗粒混合物辐射传输的影响

基于全光谱k分布(Full spectrum k distribution,FSK)模型、MIE理论和有限体积法(Finite volume method,FVM),构建了均温、均质辐射参与性气体-碳黑颗粒混合物介质热辐射传输模型,并分析了碳黑不同尺寸、不同体积浓度以及介质不同路径长度和不同温度条件下,因忽略碳黑颗粒散射所导致的介质热辐射传输特性(如辐射热流、辐射源项)的计算误差。研究结果表明:体积分数不变,增大粒径,计算误差呈现出先增大后减小的趋势;数密度不变,增大粒径,或者粒径不变,增大体积分数,均将使得计算误差相应增大;粒径、体积分数不变,增大路径长度,或者升高介质温度,均将增大计算误差。通常对于含有大颗粒、高碳黑浓度的辐射参与性气体-碳黑颗粒混合物介质,碳黑颗粒散射不能忽略。     

基于有限Fourier级数的航天器转移轨迹设计

针对航天器小推力转移轨迹的初始设计问题，利用基于三阶Fourier级数的设计方法实现了航天器小推力的多圈转移。同时，基于有限Fourier级数的形状法，对具有多个约束条件的小推力多圈转移轨迹进行了优化设计。选取了共面同轴同偏心率的初始和末端轨道位置，对所提出的方法进行了仿真验证。结果表明：与改进逆五阶多项式形状法相比，所提出的方法虽然增加了转移时间，但当转移圈数为5圈、有限Fourier级数的项数为10时，可减少将近75%的转移速度增量，同时大大减小了所需的最大推力加速度的值。     

带攻角约束的高超声速飞行器自适应反步控制器设计

针对带攻角(AOA)约束的高超声速飞行器控制问题，提出一种基于非对称时变障碍函数的非线性自适应反步控制方法。首先，将飞行器模型化为严反馈形式，以反步法为基础进行控制器设计。然后通过光滑饱和函数对名义攻角指令信号进行限幅，并保证限幅信号的可导性，限幅产生的误差通过设计辅助系统进行补偿。进而使用障碍函数对攻角指令跟踪误差进行非对称时变约束。针对不确定性和干扰，设计新型自适应律对集中干扰上界进行估计并补偿。最终通过Lyapunov理论证明了闭环系统状态量一致最终有界并且攻角始终满足时变约束。仿真结果表明，本文方法能够在满足攻角约束基础上保证良好跟踪性能。     

运载火箭推进剂复杂流动传热问题数值模拟中的模型简化方法

在运载火箭改进优化或新型运载火箭设计过程中,会遇到各种复杂的流动与传热问题。按照真实产品和物理过程直接建立的数值计算模型过于复杂,不同尺度的结构、流动与传热、长时间的飞行过程相互耦合,出现计算时间太长(数年)、收敛困难、程序调试困难等问题。针对初始条件复杂、飞行过程复杂、边界条件复杂3类具体问题,总结了3种相应的模型简化方法,分别是工程经验法、极限参数法、低维近似法,为推进剂复杂流动与传热的数值模拟提供参考。     

星载微波传感器大气温湿度物理反演算法研究

针对我国下一代静止轨道微波探测器资料应用的需求,对微波50~60GHz和183GHz的温湿遥感通道开展了不同大气条件下探测敏感高度的模拟,发现大气散射对微波探测影响明显,尤其是近地面通道,在未来数据应用时需加以注意。提出一种基于一维变分的大气温湿度廓线物理反演算法。以美国最新一代气象卫星Suomi NPP搭载的ATMS(Advanced Technology Microwave Sounder)获取的观测资料为研究对象,RTTOV(Radiative Transfer for TOVS)快速辐射传输模式为前向算子,利用一维变分算法开展了温湿度廓线反演单点试验。研究发现:地表参数对微波亮温,尤其是地面通道影响较大。同时,反演所需背景场对反演结果影响较明显,反演所得温度与背景场较接近。     

Application of improved active disturbance rejection control algorithm in tilt quad rotor

The Tilt Quad Rotor (TQR) has complex dynamics characteristics, especially in conversion mode. It is difficult to build the dynamic model of the TQR and the environmental factors have a great influence on it. To solve the problem of control in conversion mode of TQR, this paper carries out the design of the controller based on improved Active Disturbance Rejection Control (ADRC). According to the characteristics of flight in conversion mode, Tracking Differentiator (TD) with explicit model is used to solve the problem of multiple integrals when the system is high-order system. Extended State Observer (ESO) with Radial Basis Function (RBF) neural network is used to estimate and compensate for internal and external uncertainties, and the adaptive sliding mode control in Nonlinear State Error Feedback (NLSEF) is used to improve the response speed of the controller and reduce the parameters which should be tuned. Through the flight control simulation of the TQR, the validity and rationality of the control system are verified.     

Multiscale topology optimization using feature-driven method

This paper presents a multiscale design method for simultaneous topology optimization of both macrostructures and microstructures. Geometric features are extended as design primitives at both macro and micro scales and represented by Level Set Functions (LSFs). Parameters related to the locations, sizes, and orientations of macro and micro features are considered as design variables and optimized simultaneously. In the overlapping areas of different macro features, embedded microstructures are optimally figured out as the solution of the corresponding sub-optimization problem. In this study, the eXtended Finite Element Method (XFEM) is implemented for structural and sensitivity analyses with respect to design variables. This method has the advantage of using a fixed grid independent of the topology optimization process. The homogenization procedure is applied to calculate the effective properties of considered microstructures in each macro feature. Numerical examples are presented to illustrate the effectiveness of the proposed method. Results depict that the multiscale design cannot obviously improve structural stiffness compared with a solid-material design under the linear elastic condition.